NASA Technical Memorandum 112868 


//y 3 / 



Simplified Analytical Model of a 
Six-Degree-of-Freedom Large-Gap 
Magnetic Suspension System 


Nelson J. Groom 

Langley Research Center, Hampton, Virginia 


June 1997 


National Aeronautics and 
Space Administration 
Langley Research Center 
Hampton, Virginia 23681-0001 




ABSTRACT 


A simplified analytical model of a six-degree-of-freedom large-gap magnetic 
suspension system is presented. The suspended element is a cylindrical permanent magnet 
that is magnetized in a direction which is perpendicular to its axis of symmetry. The 
actuators are air core electromagnets mounted in a planar array. The analytical model 
consists of an open-loop representation of the magnetic suspension system with 
electromagnet currents as inputs. 


INTRODUCTION 

This paper develops a simplified analytical model of a six-degree-of-freedom (6DOF) 
Large-Gap Magnetic Suspension System (LGMSS). The LGMSS is a conceptual design 
for a ground-based experiment which can be used to investigate the technology issues 
associated with magnetic suspension at large gaps, such as accurate suspended element 
control and accurate sensing (ref. 1). This technology is applicable to future efforts which 
range from magnetic suspension of wind tunnel models to advanced spacecraft experiment 
isolation and pointing systems. The 6DOF model is an extension of the five degree of 
freedom (5DOF) model developed in reference 2. The suspended element is a cylindrical 
permanent magnet which is magnetized perpendicular to its axis of symmetry and the 
actuators are air core electromagnets mounted in a planar array. The electromagnet array is 
mounted horizontally with the suspended element levitated above the array by repulsive 
forces. In the nominal suspended element orientation, the axis of symmetry is horizontal 
also. The 5 DOF model developed in reference 2 was used to investigate two LQR control 
approaches for an LGMSS in reference 3. In reference 3, the simplifying assumption was 
made that the change in field and field gradients with respect to suspended element 
displacements was negligible. In reference 4 the analytical model developed in reference 2 
was linearized and extended to include the change in fields and field gradients with respect 
to suspended element displacements and the open-loop characteristics of the resulting 
system were investigated. Reference 5 developed the expanded equations (up to second 
order) for torque and force on a cylindrical permanent magnet core for two orientations of 
the core magnetization vector. One orientation was parallel to the axis of symmetry of the 
core and the other was perpendicular to this axis. In general, the higher order terms in the 
expanded equations can be neglected. However, in the case where the magnetization vector 
is perpendicular to the axis of symmetry, the expanded equations indicate that torque about 
the magnetization vector can be produced by controlling a second-order gradient term 
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directly. This allows the core to be controlled in 6DOF. In this paper the 6DOF analytical 
model is developed by following the approach detailed in references 2 and 4 using the 
equations for a cylindrical permanent magnet core uniformly magnetized perpendicular to 
its axis of symmetry which are developed in reference 5. The analytical model consists of 
an open-loop representation of the magnetic suspension system with electromagnet currents 
as inputs. 


SYMBOLS 

A system matrix (state-space representation) 

a radius of core, m 

<B input matrix (state-space representation) 

B magnetic flux density vector, T 

m matrix of field gradients, T/m 
F total force vector on suspended element, N 
F c magnetic force vector on suspended element, N 
F d disturbance force vector on suspended element, N 
F g gravitational force vector on suspended element, N 

g acceleration due to gravity (lg=9.81 m/sec 2 ), m/sec 2 

h suspension height (suspended element centroid to top plane of coils), m 
I coil current vector, A 

[IJ moment of inertia about the principal axes of the suspended element, kg-m 2 

k m constant representing magnitude of produced by in coil n 

fc ija constant representing magnitude of B ijn produced by I,^ in coil n 

*c.j)ta> constant representing magnitude of B 0j)lm produced by 1^ in coil n 

K i» *iAnax* T/A 

K* VUx> T/m/A 

^(ij)lu> ^(ij)ta/^n*x’ 

£ length of suspended element, 

M magnetization vector. A/m 

m c suspended-element mass, kg 

T total torque vector on suspended element, N-m 
T c magnetic torque vector on suspended element, N-m 
T d disturbance torque vector on suspended element, N-m 
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suspended-element rate to Euler rate transformation matrix for a 3, 2, 1 
(z, y, x respectively) rotation sequence 

inertial coordinate to suspended element coordinate vector transformation matrix 
velocity vector, m/sec 
permanent magnet core volume, m 3 
weighting matrix (eq. (36)) 
modified weighting matrix (eq. (40)) 
state vector for linearized model 

coordinates in orthogonal axis system, m 

small increment 

Euler orientation for 3, 2, 1 rotation sequence, rad 
angular velocity vector, rad/sec 

Subscripts 

electromagnet axes 
number of coils in system 
coil number 

component along x, y, z axis respectively 
partial derivative of i component in j direction 

partial derivative of ij partial derivative in k direction 
maximum value 
equilibrium condition 

Matrix Notation 

matrix 

inverse of matrix 
column vector 
row vector 

transpose of row vector 
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Dots over a symbol denote derivatives with respect to time; a bar over a symbol 
indicates that it is referenced to suspended element coordinates. 

ANALYTICAL MODEL 

This section presents a simplified analytical model of a 6DOF LGMSS which is 
developed by following the approach detailed in references 2 and 4 using the equations for 
torques and forces on a cylindrical permanent magnet core uniformly magnetized 
perpendicular to its axis of symmetry as developed in reference 5. The equations are 
simplified by using small-angle assumptions and neglecting second-order terms involving 
suspended-element motion. The permanent magnet core, or suspended element, is levitated 
over a planar array of electromagnets. Figure 1 is a schematic representation of an eight coil 
system that shows the coordinate systems and initial alignment The suspended-element 
coordinate system consists of a set of orthogonal x, y, z body-fixed axes that define the 
motion of the suspended element with respect to an orthogonal x, y, z system fixed in 
inertial space. The suspended-element coordinate system is initially aligned with the x, y, z 
system. A set of orthogonal x b -, y -, z b -axes, also fixed in inertial space, define the location 
of the electromagnet array with respect to the x, y, z system. The x b - and y b -axes are 
parallel to the x- and y-axes respectively, and the z fc - and z-axes are aligned. The centers of 
the two axis systems are separated by the distance h. The eight coil array consists of four 
coils mounted in a circular arrangement in the center and four additional coils mounted 
around the center anay. The array in the center predominantly controls the gradients of the 
fields and therefore the forces along the x-, y-, and z-axes and the outer array 
predominantly controls the magnitudes of the fields and therefore the torques about the x-, 
y-, and z-axes. The cylindrical permanent magnet core, as mentioned above, is magnetized 
perpendicular to its axis of symmetry and initially the magnetization vector is aligned with 
the positive z-axis. Therefore, control of the core involves independently controlling the x- 
and y-components of the field and their gradients in the z direction. As shown in reference 
4, independent control of B x and B n is not feasible with a single circular array of 
electromagnets (it can also be shown that independent control of B y and B^ is not 
feasible). Hence, two circular arrays with different location radii are employed. 
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Equations of Motion 


From references 2 and 4, the angular acceleration of the suspended element Q in 
suspended-element coordinates can be written as 

^ = [I C ] _1 T (1) 

pi 0 0 1 

where [l c ] = 0 I y 0 is the moment of inertia about the principal axes of the suspended 
0 0 I*_ 

element and T denotes the total torque on the suspended element A bar over a variable 
indicates that it is referenced to suspended-element coordinates. The torque T can be 
expanded as 

T = T C + T d (2) 

where T c denotes the control torque on the suspended element produced by the 
electromagnets and T d denotes external disturbance torques. The angular rates of the 
suspended element are obtained by integrating equation (1). The suspended-element Euler 
rates can be written as 

0 = [T e ]Q (3) 

where [T e ] is the suspended-element rate to Euler rate transformation matrix for a 3, 2, 1 
(z, y, x) rotation sequence. By using small-angle and rate assumptions, equation (3) 
reduces to 

6 = U (4) 

where 
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(5) 



d= 


*y 

d. 


The translational acceleration of the suspended element V in suspended-element 
coordinates can be written as 


V=(— )F (6) 

m c 

where m c is the mass of the suspended element and F denotes the total force on the 
suspended element The force F can be expanded as 

F = F C + F d + F g (7) 

where F c denotes control force on the suspended element produced by the electromagnets, 

F d denotes external disturbance forces, and F g consists of the force acting on the 

suspended element due to gravity, transformed into suspended-element coordinates. The 
suspended-element translational rates are obtained by integrating equation (6). The 

suspended-element translational rates V in inertial coordinates are given as 

V = [T m ]- 1 V (8) 

where [T m ] is the inertial coordinate to suspended-element coordinate vector- 
transformation matrix. By using small-angle and rate assumptions, equation (8) reduces to 

V = V (9) 


where 
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( 10 ) 


V = 



Magnetic Torques and Forces 


From reference 5, the torque on a permanent magnet core which is magnetized 
perpendicular to the axis of symmetry, in a given coordinate system, can be approximated 
as 


T = -vM B 

•‘cx z y 

(11) 

T cy =vM z B, 

(12) 


(13) 

where the terms that are a function of second-order gradients have been ignored for T cx and 


T cy . For simplicity define 


f a 
12 __ 4 


(14) 


Since B = [T 1 B, T and T rv in core coordinates can be written as (again using small- 

L ® J ex c j 

angle assumptions) 


T c5 = -vM ? (B y -0 2 B x +<W ( 15 > 


= vMj(B x + 0 z B y - 0 y B z ) (16) 

Obtaining T a in core coordinates is more complicated. One method of obtaining T cj is to 
transform the expanded equation for B into core coordinates using equation (A16) in 
reference 5. The gradients in core coordinates can then be calculated and substituted into 


7 



equation (40) in reference 5 and the integral over the volume taken. Following this 
approach, becomes 


v M j (c z B (xy)z + 0 x c z (B (xz)z B (xy)y ) + 6 y c z (B (zx)y B (yz)z ) 

+ ^z C z(®(yy)z — B(xx)z)) 

(17) 

From reference 5 the forces on a cylindrical permanent magnet core, for expansion of fields 
up to second order, are a function of first-order gradients only. From reference 4, the 
forces in a given coordinate system, as a function of first-order gradients, can be written as 

F c = v[<?B]M 

(18) 

where 


B xx B jy B^ 

[*]= B„ B„ 

B xz B yz B 2Z_ 

(19) 

The forces in core coordinates can be written as 


F c = v[T„][aB][T„] T M 

(20) 

For magnetization perpendicular to the axis of symmetry (along the z axis) the forces 
become 

F cX =vM f (B ]tt -0A y + 0 y (B„-B zz ) + 0 z B yz ) 

(21) 

Fey = vM I (B yz + 0 x (B a -B„) + 0^ - 0 x BJ 

(22) 

F a - = vM,(B zz -20 x B yz + 20 y B xz ) 

(23) 
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Disturbance Torques and Forces 


The assumption is made that the only significant disturbances acting on the suspended 
element is along the z - axis and is equal to its weight 


F 


g 


0 

0 

- m cS 


(24) 


where g is the acceleration of gravity. Other disturbance torques and forces are ignored. In 
suspended-element coordinates 


F.=[T„]F, (25) 

Performing the transformation (under small-angle assumptions) results in 

Fgx = 0y m cS ( 26 > 

F g y=-0x m cg ( 27 > 

Fgz = ~m c g (28) 

Linearized Equations 
The equations of motion are in the form 

( 

X = f 1 (29) 

((— )F) 
m„ J 


where 
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x T =|n, n, n s e, e, e, v„ v, v, , , zj 


(30) 


The torques and forces are functions of X and the coil currents I; thus. 


([i.r'T)] 

((— )F) 

m„ 


= f(X,I) 


where 


1 = 



(31) 


(32) 


The equations can be linearized around the nominal operating point X 0 , 1 0 by 
performing a Taylor series expansion. Neglecting second-order and higher terms and 
subtracting out X 0 results in 


where 


5X = J4SX + (B5l 


(33) 



(34) 


(35) 


and 
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[j 0 0 0000 0 0 000 

1/I 7 0 0000 0 0 000 

0 1/Ij 0000 0 0 000 

0 0 1000 0 0 000 

0 0 0100 0 0 000 

0 0 0010 0 0 000 

0 0 0 0 0 l/m c 0 0 0 0 0 

0 0 0 0 0 0 l/m c 0 0 0 0 

0 0 0 0 0 0 0 l/m c 0 0 0 

0 0 0000 0 0 100 

0 0 0000 0 0 010 

0 0 0000 0 0 001 

Expanding J? results in 

'dTj/dQj ... &T x /dz 

^I 7 /dQ- % <?T 7 /<?Q 7 ... 
y? = W, ... (37) 

dv-jdQj ... dv^/dz 

which reduces to 
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Finally, by using the expressions for torques and forces developed earlier (eqns. (15)-(17), 
(21)-(23), and (26)-(28)), A becomes 
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vM z / 1 5 0 0000 0 0 0 000 

0 vMj/I y 0000 0 0 0 000 

0 0 vM z /Ij 0 0 0 0 0 0 0 0 0 

0 0 0100 0 0 0000 

0 0 0010 0 0 0000 

0 0 0001 0 0 0000 

0 0 0 0 0 0 vMj / m c 0 0 000 

0 0 0 0 0 0 0 vMj / m c 0 0 0 0 

0 0 00000 0 vM f /ra c 0 0 0 

0 0 0000 0 0 0100 

0 0 0000 0 0 0010 

0 0 0000 0 0 0001 


Next, expanding <2 results in 

'3 r x ldi x #T t /d I 2 ... dTj/dL m ' 

#r y /d i, &r y /di 2 ... 

<B = W 1 ... (41) 

Evaluating the first term in equation (41) at X,, results in 

(42) 

Since the fields and gradients are linear functions of coil currents, the components of B y 
produced by coil n of an m-coil system can be written as 

(43) 

where is the maximum coil current, is a constant that represents the magnitude of 

B yn produced by 1^, and 1^ is the coil current. To simplify, define 
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( 44 ) 


K yn = ^max 


For the total system, B can be written as 


B y =|K,Jl 


(45) 


where 


K,=|K„ K* - V| 


(46) 


and I is defined by equation (32). Since the elements of ]_K y J are constants 


<7B y / <21 = |_K y J 


(47) 


Similar results are obtained for the other fields and gradients. Terms in B related to the 
identities 6 X - 6 y - x = V-, y = V y , and z = V- are zero. Then © 

becomes 


© = W, 


"L K J 

I KxJ 

'zl_K(iy)2 J 

Loj 

Loj 

Loj 

L*J 

W 

L K *d 

Loj 

LOJ 

LOJ 


(48) 
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Initial Conditions 


The suspended element is assumed to be initially suspended in equilibrium at a distance 
h above the electromagnet array with the suspended-element coordinates initially aligned 
with the inertial coordinates as shown in figure 1. In equilibrium, F x = F y = 0 and 


F*=m c g 

From equations (15)-(17) and (21)-(23), we have 

®x = = ®(xy)z = ®xz = ®yz = ® 


(49) 


(50) 


and 


B 22 =^- (51) 

vMj 

In equilibrium, by using the relationship of equation (51), elements (7, 5) and (8,4) of the 
matrix in equation (39) reduce to 


m c g 

vM 5 


+ (B„ — ) = B xx 


(52) 


and 


^-V-^r = - B yy 


vM x 


(53) 


From equation (45), the controlled fields and gradients as a function of I,, can be written as 
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( 54 ) 


"B x ' 


IKJ " 



L k J 

B« 


L K «J 



L k »J 

K 


L K J 

_B(xy) Z _ 


|K, W J 


I,, can be found by inverting the K matrix in equation (54) using the generalized inverse. 
This produces a solution where the 2-norm of the current vector is minimized (i.e., 

minimum £i 2 ). Once I c is determined, the uncontrolled fields and gradients required to 
complete the A matrix can be calculated. 


As noted in references 3 and 4, one of the objectives of the LGMSS development is to 
allow positioning of the suspended element through large angles in yaw (0 Z ) up to 360°. As 

the suspended element is rotated, the equilibrium currents will change. Following the 
approach detailed in the Appendix of reference 2, the equilibrium currents can be developed 
as a function of yaw angle and initial torques and forces on the suspended element 

Assuming only yaw displacement [T m ] becomes 


[T„] = 


"c 6 Z s 6 Z 
-sd z c6 z 
0 0 


O' 

0 ! 

1 


(55) 


Since B = [T m ]B, from equations (11) and (12) 


(s0 z B x -c0 z B y )1 

.(c0 z B x + s0 z B y )J 

Using the approach discussed earlier to obtain equation (17), T cz becomes 
T C z = y M z c z ((c 6 Z — s 0 z )B (xy)z + c0 z s0 z (B (yy)z — B (xx)z )) 



(56) 


(57) 
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Simplifying equation (57) results in 


T a = vM 1 c 1 (c20,B w) ,+|s20 1 (B,„„ -B ( „»)) (58) 

The forces, from equation (20) become 

’ F c*l [’(c0 2 B ][2 +s0 z B y2 )~ 

F c y =VM, (ce.Bj.-^A.) (59) 

[ B 2z 

In terms of yaw angle and coil currents, the torques and forces become 

(se,|Kj-c«,|K,J) 

(ce,LK,J+sfl,L K .J> 

c,(c2e,[K w „J + i S 2e I ([K,„ )1 J - |k w * J» 

(c«,L K =J+ se ,L K ^J> 

(c9,L k pJ -s9 >L k =J) 

|kJ 

CONCLUDING REMARKS 

This paper has developed a simplified analytical model of a six-degree-of-freedom large- 
gap magnetic suspension system. The suspended element is a cylindrical permanent magnet 
that is magnetized perpendicular to its axis of symmetry and the actuators are air-core 
electromagnets mounted in a planar array. The analytical model is an open-loop 
representation with electromagnet currents as inputs. The model should be useful in 
analyses and simulations, in the development of control system approaches, and in 
evaluations of overall system performance. 
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Figure 1 . Initial alignment of coordinate systems for 6 DOF 
Large-Gap Magnetic Suspension System. 
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